The Breakdown of Kinetic Theory in Granular Shear Flows 
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We examine two basic assumptions of kinetic theory- binary collisions and molecular chaos- using 
numerical simulations of sheared granular materials. We investigate a wide range of densities and 
restitution coefficients and demonstrate that kinetic theory breaks down at large density and small 
restitution coefficients. In the regimes where kinetic theory fails, there is an associated emergence 
of clusters of spatially correlated grains. 
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For granular materials, kinetic theory has been the 
primary strategy used to systematically derive hydrody- 
namics equations, starting from elementary assumptions 
about grain-grain interactions Q. This has led to much 
interest in applying predictions from the theory to real- 
istic granular flows (for reviews, see Q), and recent work 
continues in this direction However, kinetic theory 

strictly applies only to dilute gases, and the extent that 
it applies to the dense regime remains unclear. 

In this Letter we perform tests of the fundamental as- 
sumptions of kinetic theory, using the Contact Dynamics 
(CD) algorithm. We find that kinetic theory is severely 
limited by the assumption that only binary interactions 
occur between grains. Instead, an effective many body in- 
teraction arises that is a direct consequence of persistent 
contacts in the dense regime. In Fig. ^ we characterize 
the failure of the binary collision assumption, using spa- 
tial force correlations to approximate the average number 
of grains N c that form a cluster in contact. As we see, 
the cluster size increases when going to low restitution 
coefficient and high density, which should limit the rele- 
vance of kinetic theory. This Letter provides quantitative 
estimates of this breakdown. 

Most kinetic theory research starts with the Boltzmann 
equation, which is derived from the BBGKY hierarchy, 
and then finds its solutions 0, 0, 0. However, certain 
assumptions are necessary to derive the Boltzmann equa- 
tion. We begin by discussing two of these assumptions. 

Consider a system of N grains and the evolution equa- 
tion for the iV-body probability distribution function 
(pdf) f( N >(r, p), where (r, p) denotes the set of all po- 
sitions and momenta for the system (with the notation 
r = {ri}). This equation is simply a statement of con- 
servation of probability and reads 



d f(N) 



dt 
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(1) 



where we have decomposed the force on each grain as a 
sum over pairs: i*j = J2j 

The BBGKY hierarchy is derived from Eq. (JTJ by inte- 
gration 6J. This hierarchy is the set of N— 1 equations for 
the evolution of the n-body pdfs f^ n \ with 1 < n < N. If 
the force resulting from a pair interaction depends only 



on the positions and velocities of the interacting pair, 
then the evolution equation of each depends only on 
j(n+i) This ig thc classical form of thc BBGKY hierar- 
chy. 
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FIG. f : Contours of the average cluster size N c as a function 
of restitution coefficient e and packing fraction v. Kinetic 
theory applies to binary interactions between grains (N c = 2) . 

When applying the derivation of the BBGKY hierar- 
chy to granular materials, additional care must be taken. 
For dilute hard-sphere gases the derivation applies since 
interactions result from binary collisions and thus pair 
forces depend only on the positions and velocities of the 
interacting pair. However, in dense systems of dissipative 
grains the situation is different. In this case there may 
be clusters of grains in persistent contact, as illustrated 
in Fig. [2 If clusters have formed, the force between any 
pair in the cluster will depend not only on the positions 
and velocities of the interacting pair, but also on the po- 
sitions and velocities of all other grains in the cluster. In 
this case, the hierarchical structure of the BBGKY equa- 
tions is not guaranteed. Therefore, in order to derive 
the BBGKY hierarchy for granular materials, we must 
make the binary collision assumption. This stipulates 
that only binary collisions occur, thereby assuring that 
pair forces depend only on the positions and velocities of 
the interacting pair. 

The Boltzmann equation follows from the first equa- 
tion (n = 1) of the BBGKY hierarchy, which relates 
to f( 2 K A second assumption is also required, the molec- 
ular chaos assumption, which simplifies this equation by 
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setting 

.f {2) (r 1 ,r 2 ,vi,v 2 ) = x( r i, r 2)f W {ri,vi)f {1) (r2,v 2 ), (2) 

where x describes possible correlations in the positions of 
particles. This assumption allows us to approximate the 
first equation of the BBGKY hierarchy as a non-linear 
equation for /W: this is the Boltzmann equation. 

The Boltzmann equation relies on the assumptions of 
binary collisions and molecular chaos. In the rest of this 
Letter we test these fundamental assumptions, using two- 
dimensional CD simulations 7] of frictionless granular 
materials in simple shear flow at constant volume. The 
simulations are performed using Lees-Edwards boundary 
conditions, which ensure translational invariance. The 
density, restitution coefficient, and shear rate are pre- 
scribed and other observables are measured. 

In Fig. we show two representative screenshots from 
our simulations, in steady state, for identical shear rate 
and density, but different restitution coefficients, e = 0.92 
and e = 0. A small time interval is chosen, and in both 
cases grains that collide during this time interval are col- 
ored. Different colors corresponding to separate contact 
networks. For e = 0.92 the interacting grains are well 
spaced and tend to occur in pairs, whereas for e = 
the interacting grains tend to form large clusters. These 
clusters indicate the emergence of persistent contacts for 
small values of e (high inelasticity). 




FIG. 2: (Color online). Steady state screenshots of sheared 
granular materials at two different restitution coefficients: e = 
0.92 (left) and e = (right). Grains involved in a collision 
during a small time period are colored, with different colors 
denoting different contact networks. Interactions are binary 
for e = 0.92 whereas large clusters form for e = 0. 

We now test the binary collision assumption, which is 
primarily a statement about the forces between pairs of 
grains. The idea is to measure the relative contribution 
to momentum transport from binary collisions and from 
clusters in persistent contact. A quantitative measure of 
these contributions can be obtained by measures of the 
stress tensor. The "static" stress X s reads: 

= \ £(A + £>i)fiUaF U/ 5 (3) 

where a, (3 denote components and i, j denote grains, 
Di is the diameter of grain i, ny the unit normal vector 



at contact between the pair and A is the area of 

the simulation cell. This quantity measures the true mo- 
mentum transport in a microcanonical configuration. At 
each time step, the CD algorithm determines the contact 
forces Fy by upholding constraints relevant to perfectly 
rigid contact between grains @- 

When a binary collision occurs, the final relative veloc- 
ity of the colliding pair is set equal to the initial relative 
velocity, multiplied by — e. The CD algorithm calculates 
the force between the pair based on the instantaneous im- 
pulse that produces this final relative velocity. Because 
of the time-discretization, this "binary collision force" 
is approximated over a time interval At by a constant 
force equal to the instantaneous impulse divided by At: 
Fy c = ^| (1 + e) (vi — vj) • ny ny, where fj^j is the re- 
duced mass of grains i and j, and Vj is the pre-collisional 
velocity of grain i. However, when multiparticle collisions 
occur, the total forces Fy differ from the binary collision 
forces Fy hc . Replacing Fy by Fy hc in Eq. © provides 
the flux of momentum that would be transported if all 
forces resulted from binary collisions: 

^%A = -2^X)pii(A+i5i)fi«ofi«|9(vi-vj)-fiy. (4) 

We call this tensor the "collisional" stress tensor: it is 
defined at any time, even in the presence of multi-contact 
interactions, and is an approximation to the static stress. 

A theory that assumes binary collisions and is capable 
of taking into account all correlations and providing an 
exact expression for the distribution of velocities between 
incoming pairs of grains would only account for Xjjjj, 
but never X* g. Because most kinetic theories assume 
binary collisions, the core question is whether X^g is a 
reasonable approximation to X* g . 

To answer this question, we further decompose stresses 
into pressure p and shear stress s. Pressure is one-half of 
the trace of the tensor and the shear stress is defined as 
either of the off-diagonal elements of the symmetric stress 
tensor. In Fig.|21we plot data from our simulations for the 
static pressure divided by the collisional pressure p s /p bc 
and the static shear stress divided by the collisional shear 
stress s s /s bc as a function of packing fraction, for a vari- 
ety of restitution coefficients. For restitution coefficients 
near unity and relatively low packing fraction, the static 
values are equal to the collisional values and the ratios in 
Fig. [3] are close to unity. However, for large packing frac- 
tions and small restitution coefficients, the static values 
become larger than the collisional values. This signals 
the breakdown of the binary collision assumption. 

This first numerical test quantitatively demonstrates 
that the collisional stress tensor is not an adequate ap- 
proximation of the true static stress tensor in certain 
regimes of granular shear flow. It rules out the possibility 
that a theory based on the binary collision assumption 
can be applied to predicting the static stress at high den- 
sity or low restitution. 
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FIG. 3: The static pressure divided by the collisional pressure 
(p s /p bc ) and the static shear stress divided by the collisional 
shear stress (s s / s hc ) as a function of packing fraction, for a 
variety of restitution coefficients e. Values larger than one 
represent the breakdown of the binary collision assumption. 

Next, we examine the molecular chaos assumption: we 
test whether a broad array of kinetic theories succeed 
in accounting for the collisional stress. In two dimen- 
sions, kinetic theories that assume both binary collisions 
and molecular chaos make a prediction for the collisional 
pressure p kt 

P kt = (1 + e)x»V* (5) 

where v is the packing fraction and \ is the pair correla- 
tion function at contact. This prediction is proportional 
to p* — nm5v 2 /2 where n is the number density, m is 
the average mass, and Sv 2 is the average square of the 
fluctuating velocity (the granular temperature). 

Our second numerical test compares the prediction for 
the collisional pressure to the actual collisional pressure 
measured in the simulations. We determine all param- 
eters in Eq. (JSJ directly from the simulations: e and v 
are prescribed, x ancl P* are measured. Following other 
studies |I| we measure x using the collision frequency u> 
and the following formula from kinetic theory: 

u> = V2Tr5v 2 ~x n<J (6) 

where a is the average grain diameter. 

Using Eqs. © and © we measure, without any fitting 
parameters, the approximation to the collisional pressure 
resulting from the molecular chaos assumption. This is 
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FIG. 4: Main figures: The static pressure p s , collisional pres- 
sure p bc , and kinetic theory prediction p kt , all normalized by 
p* , for e = 0.92 (top) and e = (bottom). Insets: Velocity 
correlations at contact, normalized by the granular temper- 
ature Sv 2 . p kt from Eq. JS) overestimates p bc when velocity 
correlations are large and positive. 

reported in solid lines on Fig. where we have also plot- 
ted raw data for both the collisional and static pressure. 
For e = 0.92 there is excellent agreement between the ki- 
netic theory prediction and the collisional pressure, even 
for large values of packing fraction. For e — the molecu- 
lar chaos assumption leads to an overestimate of pressure 
at all packing fractions. 

We expect this overestimate to result from correlations 
of the pre-collisional velocities: if the velocities of two in- 
coming grains are positively correlated then their relative 
velocity is smaller, and the collisional pressure is thereby 
reduced. Because the molecular chaos assumption does 
not incorporate these correlations, it overestimates the 
collisional pressure. 

The insets of Fig. 0] contain measurements of the pre- 
collisional velocity correlations (5vi6v2), where Sv is the 
fluctuating part of the velocity in the direction paral- 
lel to the vector connecting the grain centers ny. Wc 
measure correlations in this direction because only these 
velocity components contribute to the collisional stress in 
Eq. (J2J. The average is performed over a disk centered 
on one grain with a radius of 1.8cr, although the results 
do not depend on the size of the averaging disk. Pairs of 
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FIG. 5: Inset: Plots of the spatial force- force correlations 
|C(£)| (in arbitrary units on a log scale) for e=0 and three 
packing fractions v — 0.7 (bottom) v = 0.77 (middle) and 
v — 0.81 (top). Main Figure: The average cluster size N c as a 
function of packing fraction, for many values of the restitution 
coefficient e. 

grains that collided in the previous time step are excluded 
from the average in order to ensure the correlations are 
truly pre-collisional. We observe that large values of the 
correlation correspond to packing fractions where kinetic 
theories based on the molecular chaos assumption overes- 
timate the collisional pressure. However, the value of the 
correlation does not correspond to whether the collisional 
stress is a good approximation to the static stress- this 
seems to be related only to the breakdown of the binary 
collision assumption. 

In order to further understand how the binary collision 
assumption breaks down, we return to the observation 
that dense granular flows organize into clusters of inter- 
acting grains, as illustrated in Fig. The failure of the 
binary collision assumption is related to the formation of 
these clusters and the fact that the stress tensor is not 
determined solely by two particle interactions. 

A cluster of grains may be defined as a region in the 
material over which forces are correlated. To deter- 
mine the average number of grains N c in these correlated 
clusters, we measure the spatial force-force correlations 
C(£) = (F(0) ■ F(£)), where £ is a positive distance mea- 
sured in grain diameters and F is the total vector force 
acting on a grain. If a grain is isolated, so that there 
are no forces acting on it, it is not included in this aver- 
age. We then define N c proportional to the square of the 
correlation length 

/ £C(£)d£ 



(*> = 



JC(£)d£ 



(7) 



and normalize so N c = 2 for e = 0.92 and low density. 
We choose this normalization because we have observed 
(see Fig. [3J| that the binary collision assumption, which 
corresponds to N c =2, is appropriate for dilute, nearly 
clastic granular materials. 

Our measurements of N c are presented in Fig. \5\ along 
with measurements of C(£). The force correlations fluc- 



tuate at small distance and exhibit an exponential de- 
cay at large distance. The values of {£) determined 
from Eq. (JJJ match the exponential decay (when we plot 
e -^/w) f or the densities and restitution coefficients we 
have investigated. We notice from Fig. [S] that the diver- 
gence of N c close to jamming nicely echoes the divergence 
of pressure and shear stress ratios in Fig. [3] This con- 
firms that the formation of clusters is directly related to 
the breakdown of the binary collision assumption. 

The measurement of N c allows us to partition the 
phase space of granular shear flow into regions where ki- 
netic theory applies (N c = 2) and regions where it does 
not (N c > 2). In Fig. ^ wc plot contours of N c as a 
function of restitution coefficient and packing fraction. 
Although numerical noise prevents us from plotting the 
contour N c = 2, Fig.^provides an estimate of the regime 
where kinetic theory applies. 

We have presented two numerical tests of the funda- 
mental assumptions of kinetic theory in granular materi- 
als: first we have observed the breakdown of the binary 
collision assumption for large densities and small restitu- 
tion coefficients; second we have demonstrated that the 
molecular chaos assumption is not valid for small resti- 
tution coefficients, due to pre-collisional velocity correla- 
tions. In order for an approach based on kinetic theory 
to be useful at high density, the deficiencies in these core 
assumptions must be addressed. Although the molecular 
chaos assumption can in principle be addressed by incor- 
porating velocity dependent terms in Eq. @, it seems to 
us that the failure of the binary collision assumption may 
be much more difficult to overcome in a standard kinetic 
theory. Successful theories of granular materials in the 
dense regime must incorporate clustering. 
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